Novel schemes for measurement-based quantum computation 
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We establish a framework which allows one to construct novel schemes for measurement-based quantum 
computation. The technique further develops tools from many-body physics - based on finitely correlated or 
projected entangled pair states - to go beyond the cluster-state based one-way computer. We identify resource 
states that are radically different from the cluster state, in that they exhibit non-vanishing correlation functions, 
can partly be prepared using gates with non-maximal entangling power, or have very different local entangle- 
ment properties. In the computational models, the randomness is compensated in a different manner. It is shown 
that there exist resource states which are locally arbitrarily close to a pure state. Finally, we comment on the 
possibility of tailoring computational models to specific physical systems as, e.g. cold atoms in optical lattices. 

PACS numbers: 03.67.-a, 03.67.Mn, 03.67.Lx, 24.10.Cn 



No classical method is known which is capable of effi- 
ciently simulating the results of measurements on a general 
many-body quantum system: the exponentially large state 
space renders this a tremendously difficult task. What is 
a burden to computational physics can be made a virtue in 
quantum information science: It has been shown that multi- 
particle quantum states can form resources for quantum com- 
puting Indeed, universal quantum computation is possible 
by first preparing a certain multi-partite entangled resource - 
called a cluster state Q], which does not depend on the algo- 
rithm to be implemented - followed by local measurements 
on the constituents. This idea of a measurement-based "one- 
way computer" (QCc) lUt] has attracted considerable attention 
in recent years. Progress has indeed been made concerning a 
systematic understanding of the computational model of the 
one-way computer as such fllj, |g, 0> [Ml ■ Quite surprisingly, 
this contrasts with the lack of development of new compu- 
tational models or novel resource states beyond that original 
framework. To our knowledge - no single model distinct from 
the QCc has been developed based on local measurements on 
a fixed, algorithm-independent qubit resource state. Hence, 
questions of salient interest seem to be: Can we systemati- 
cally find alternative schemes for measurement-based quan- 
tum computation? What are the properties that distinguish 
computationally universal resource states? 

These questions are clearly central when thinking of tai- 
loring resource states to specific physical systems, e.g. to 
cold atoms in optical lattices, purely linear optical systems or 
condensed-matter ground states. They are also of key interest 
when addressing the question what flexibility one has in the 
construction of such schemes, and what properties of may ulti- 
mately be relaxed. The problem is also relevant to many-body 
physics, when the question of efficient classical simulatability 
121 is addressed: Quantum states may thought of being or- 
dered according to their computational potency, universal and 
efficiently simulatable states forming the respective extremes. 

In this work, we demonstrate how methods from many- 
body physics can be extended to develop schemes for 
measurement-based quantum computation (MBC). Starting 
from the concepts of matrix- prod uct, projected entangled pair, 
and finitely correlated states lldllTlll . we develop a framework 



broad enough to allow for the construction of novel univer- 
sal resources and models. The notion of universality in the 
context of one-way computing was recently addressed in Ref. 
13. A universal resource in their sense is a family of states 
out of which any other state can be obtained by local measure- 
ments on a subset of sites. It follows from the definition that 
many states cannot be universal: E.g. states which are locally 
non-maximally entangled, have non-vanishing two point cor- 
relation functions (OjOj+r) — (Oi) (Oi+ r ) or a non-maximal 
localizable entanglement 11411 . Complementary to this ap- 
proach, we refer to a device as a universal quantum computer, 
if it can efficiently predict the outcome of any quantum algo- 
rithm. A state will hence be called a universal resource if one 
can, assisted by the results of local measurements on the state, 
efficiently predict the result of any quantum computation. 
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FIG. 1: Two resources for universal measurement-based quantum 
computing. Fig. la depicts a weighted graph state, where solid lines 
correspond to a controlled 7r-phase gate, dashed lines to n/2. Fig. 
lb represents a scheme deriving from an AKLT-type model. Dashed 
lines represent a state with non-vanishing correlation functions, solid 
lines correspond to 7r-phase gates in the {|1), |2)}-subspace. 

To exemplify the power of our framework, we describe 
three new models for MBC in quantum mechanical lattice sys- 
tems. In all these models, the randomness is compensated in 
a manner different from the QCc- They are chosen to high- 
light that, intriguingly, many properties of the original one- 
way computer may be relaxed, while still retaining a universal 
model for quantum computing: (i) We find resources exhibit- 
ing non-vanishing two-point correlations (which are typical 
for natural many-body ground states). The original discussion 
of the QCc depended on the fact that the cluster can be pre- 
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pared by mutually commuting unitaries (technically a quan- 
tum cellular automaton (QCA) |fl2ll ). Commutativity enables 
one to logically break down a QCc -calculation into small 
parts corresponding to individual gates; however, it implies 
severe restrictions, such as that the correlations vanish com- 
pletely outside some neighborhood. Hence our framework can 
prove universality for states not amenable to any QCA-based 
technique, (ii) We treat a universal weighted graph state with 
partly weakly entangled bonds, (iii) In the final part - using 
different techniques - we present a family of states which are 
universal, yet are locally arbitrarily close to a pure state. 

Matrix product states. - Starting point in the 1-D setting 
is the familiar notion of a matrix product state (MPS) fioll . 
We will first look at the simple case of a chain of n qubits. 
Its state is specified by (i) an auxiliary D-dimensional Hilbert 
space, called correlation space, (ii) two operators A[0],^4[l] 
on C D , and (iii) two vectors \L), \R) representing boundary 
conditions. One has explicitly 



|tt)= (R\A[s n ]...A[s 1 }\L)\s 1 ,. 
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In order to generalize Eq. (Q3 to 2-D lattices, we need to 
cast it into the form of a tensor network. Setting Li = 

(i\L),A[s]ij = (j\A\i), we arrive at (s 1 , . . . , s„|*) = 
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Computational tensor networks. - While the 1-D setting is 
awkward enough, the 2-D equivalent is completely unintelli- 
gible. To cure this problem, we introduce a graphical notation 
111511 which enables an intuitive understanding beyond the 1 -D 
case. Tensors will be represented by boxes, indices by edges: 



L r =\L\ — , A[s)i 



A[s] 



Rh 



A single-index tensor can be interpreted as the expansion co- 
efficients of either a "ket" or a "bra". Sometimes, we will in- 
dicate what interpretation we have in mind by placing arrows 
on the edges: outgoing arrows designating "kets", incoming 
arrows "bras". Connected arrows designate contractions of 
the respective indices. If \<p) is a general state vector in C 2 , 
we abbreviate (0|O) A[ti\ + <0| 1> A[l] by A[<j>}. The overlap 
of with a local projection operator is easily derived: 



<&l)l*> = 



L 




A[<t>x] 




A[K] 







Eq. d2J should be read as follows: Initially, the correlation 
system is in the state \L). Subsequent measurements of local 
observables with eigenvectors \<j>i) at the i-th site induce the 
evolution A[<^], thereby "processing" the state in the correla- 
tion space. The probability of a certain sequence of measure- 
ments to occur is given by the overlap of the resulting state 
vector with \R). An appealing perspective on MBC suggests 
itself: Measurement-based computing takes place in corre- 
lation space; the gates acting on the correlation systems are 
determined by local measurements. The crucial new insight 
compared to previous treatments of MPS and PEPS in the 
context of many-body physics [II]] or MBC JH] is that 



the matrices used in the parametrization of an MPS can be 
directly understood as quantum gates on a logical space. We 
will refer to this representation of MBC, as a computational 
tensor network (CTN). 

The graphical notation greatly facilitates the passage to 2-D 
lattices. Here, the tensors A[0/1] have four indices A[s]i : d, r ,u, 
which will be contracted with the indices of the left, right, 
upper and lower neighboring tensors respectively: 
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for various boundary conditions L, D, R, U. Notably, it is 
known [11] that classical computers cannot efficiently per- 
form the contraction appearing in Eq. (0). This fact is an 
incarnation of the power of quantum computers and no prob- 
lem to our approach. We will now describe several examples, 
demonstrating the versatility of our framework and showing 
how - surprisingly - many reasonable assumptions about uni- 
versal resources turn out to be unnecessary. In what follows, 
we use the standard notation X, Y, Z for the Pauli operators, 
H for the Hadamard gate and S — diag(l,i) for the tt/4- 
gate. The controlled <p-phase gate is |0, 0)(0, 0| + |0, 1) (0, 1| + 
|1, 0)(1, 0| + e*|l, 1|. Lastly, |±> = 2- 1 / 2 (|0) ± |1}). 

Two-point correlations. - Here we consider a resource with 
exponentially decaying correlation functions, in a way as it 
occurs in ground states, but not in states resulting from a QCA. 
To be brief, we first describe a 1-D setting, turning to 2-D 
structures later. Define G := exp(iir/kX) for some integer 
k > 2. The relations 



A[s 



= G\s) r ( S \i, \L) = \+), \R) = G-i\+) 



define a state vector \ for a chain of qubits. The two-point 
correlations never vanish completely: One finds ifioll that 

(ZiZi+n+l) — (Zi) (Zi +n+ i) oc 

where £ = 2sin 2 (7r//c) — 1. Still, all single-qubit unitaries 
on the correlation system can be realized by local physical 
measurements. Ignoring global factors (as we will do when 
possible), one computes: 



A\ 



G, 



GZ; 



A[X] 



GZ X , 



where the r.h.s. is a compact notation for the two equations 
on the left: An observable as the argument to A[] denotes a 
measurement in the corresponding eigenbasis. The outcome 
of the measurement is assigned to a variable; here x = in 
case of the +l-eigenvalue and x = 1 in case of —1. Local 
X-basis measurements hence cause the state of the correla- 
tion system to be transported from left to right (up to local 
unitaries). When measuring several consecutive sites in the 
X-basis, the overall operator applied to the correlation system 
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is given by B := . . . GZ X2 GZ Xl . Assuming that we intended 
to just transport the information faithfully, we conceive B as 
an unwanted by-product. To understand this structure, con- 
sider the following elementary statement: Let A,B be ma- 
trices having finite order flUl . Every element in the group 
generated by A, B can be written as AB kl AB k2 A . . . AB kn 
for some n G N and fcj 6 {0,1}. Applied to our situa- 
tion: The by-product operators form a finite group B gen- 
erated by G, Z. The grou p property gives a possibility to 
cope with by-products 11911 : Assume that at some point the 
state vector of the correlation system is given by B\ijj), for 
some unwanted B G B. Transferring the state along the 
chain will introduce any by-product B' G B after a finite ex- 
pected number of steps. In particular, B' = B^ 1 will oc- 
cur, leaving us with B~ 1 B\ip) = \ip). Note that this tech- 
nique is completely general: it can deal with any finite by- 
product group (see further examples below). Moving on, a 
measurement in the 2~ 1 / 2 (|0) ± e 4 *|l) basis induces the op- 
erator S(4>) := diag(l, e 1 ^) on the correlation system (up 
to by-products). But by the preceding discussion, we can 
also implement GS{4>)G~ 1 , which is sufficient to generate 
any single-qubit unitary JTt]. Lastly, it is easy to see that Z- 
measurements prepare a known state in the correlation system 
and conversely can be used to read it out. 

AKLT-type states. - In this example, we consider ground 
states of nearest-neighbor spin-1 Hamiltonians of the AKLT- 
type, as they are well-known in the context of condensed- 
matter physics llOll . We investigate ground-states induced 

-1/2 



by ^[0] 



A[2] 



*■ = H, 
= (X + iY)2 



A[l] 



-1/2 



= (X - iY)2' 

This is the exact unique 
ground state of a nearest-neighbor frustration-free gapped 
Hamiltonian (in the original AKLT model H is replaced by 
Z in the definition of A[0]). Proving that any single-qubit uni- 
tary can be realized on the correlation space commences in a 
similar way as before. Measurements in the {|0}, 2~ 1 ' 2 (|1) ± 
e^|2))}-basis gives rise to H, S((f>) or ZS{(j>), depending on 
the measurement outcome. The finite by-product group is in 
this case generated by H, Z. But that is all we need to show, 
as gates of the form S(<j)), HS(4>)H generate all of SU{2). 

Weighted graph states. - Both previous examples can be 
embedded into 2-D lattices, universal for computation (see 
Fig. lb). A general technique for coupling 1-D chains to 2- 
D universal resources will be discussed by means of a further 
example: the weighted graph state JH O shown in Fig. 1 
(a). In the figure, vertices denote physical systems initially 
in |+), solid edges the application of a controlled 7r-phase 
gate and dashed edges controlled 7r/2-phases, so some of 
the entangling gates do not have maximal entangling power. 
The resource's tensor representation (acting on a D = 2- 
dimensional correlation space) is given by 



A[s 



= S S \+) ru S s \+) lu Z s \+) r (s\lMrMl, (4) 



where s G {0, 1}. Indices are labeled ru for "right-up" to Id 
for "left-down". Boundary conditions are |0) for the ru, lu, re- 
directions; |+) otherwise. The broad setting for our scheme 



is the following: the correlation system of every second hor- 
izontal line in the lattice is interpreted as a logical qubit. In- 
termediate lines will either be measured in the Z-eigenbasis 
- causing the logical bits to be isolated - or in the y -basis - 
mediating an interaction between adjacent logical qubits. 

We will first describe how to realize isolated evolutions of 
logical qubits. According to Eq. (O the tensors A[0/1] fac- 
tor, allowing us to draw only the arrows corresponding to the 
factors of interest; so e.g. A[s] ->■ = Z s \+) r . We find 



A[Zi_i, u ] 








A[Xi] 












A[Zi- XA ] 











HS ZXi+Zi , (5) 



where = 2f_i )tt + z%-\,d + z l+1 . u + Zi+i.d- Eq. (0 is 
of the kind treated before in the case of 1-D chains. In- 
deed, using the same techniques, one sees easily that general 
local unitaries can be implemented by measurements in the 
2- 1 / 2 (|0)±e^|l)) basis. The by-product group here is given 
by the full single-qubit Clifford group. Turning to two-qubit 
interactions, consider the schematics for a controlled-Z gate 
(we suppress adjacent sites measured in the Z-basis): 



A[X] 



A[X] 



A[X] 



A\Y] 



A[X] 



A[X] 



(6) 



A[X] 



In detail: first we perform the X-measurements on the sites 
shown and the Z-measurements on the adjacent ones. If any 
of these measurements yields the result "1", we apply a Z- 
measurement to the central site labeled Y and restart the pro- 
cedure three sites to the right lfl9tl . If all outcomes are "0", a 
y-measurement is performed on the central site, obtaining the 
result y. Let us analyze the gate step by step. For c G {0, 1}, 



A[X] -*] A[X) [*\ A[X] 
% — 



HZ C 



S c 



^ 



A[Y] 



( 1 + { _ ir +y iS c®S°)\+) lu \+) ru 



S c \+) S c \+) 



|c) J A[X] U A[X] U A[X] U ~ SC \ + ^u SC \+)m H \ c ), 



In summary, the evolution afforded on the upper line is H ( 1 + 
(-l) c+v iZ) oc HSZ y+c , equivalent to Z c up to by-products. 
This completes the proof of universality. For completeness, 
note that we never need the by-products to vanish for all log- 
ical qubits simultaneously. Hence the expected number of 
steps for the realization of one- or two-qubit gates is a con- 
stant in the number of total logical qubits. 

Entanglement properties of universal resources. - In this 
section, we further investigate - using different methods - to 
what extent the entanglement properties of the cluster state 
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can be relaxed while retaining universality in the above sense. 
More specifically, we ask whether one can find qubit resources 
that are (i) universal for QCc, (ii) translationally invariant, 
(iii) which have an arbitrarily small local entropy and localiz- 
able entanglement (LE), and (iv) from which not even a Bell 
pair can be deterministically distilled? 

To show that - rather surprisingly - this is indeed the case, 
we will encode each logical qubit in blocks of 2k + 1 horizon- 
tally adjacent physical qubits. Here, k is an arbitrary param- 
eter. The first k qubits per block will take the role of "code- 
words", the final k + 1 are "marker qubits" used in a construc- 
tion to make the resource translationally invariant. We start by 
preparing a regular cluster state in the respective first qubit of 
each block. Then, we encode the states of each of these first 
qubits according to |0) h-> \O k ) := |0)® fe and |1) h-> \W k ) := 
fc-V 2 (|0,...,0,l) + |0,...,1,0>+ |1,...,0,0». The rear 
k + 1 qubits of each block are prepared in |0, . . . , 0, 1). Call 
the resulting state vector \<fi). Finally, the translation invari- 
ant resource we will consider is >J/) = Ylt=o ^l^)' where T 
implements a cyclic translation of the lattice in the horizontal 
direction. To realize universal computing, we pick one block 
and measure each of its qubits in the Z-basis. In this way, one 
can deterministically distinguish the states T'|^>) correspond- 
ing to different values of t. Indeed, the maker codewords 
guarantee that a sequence of k "0"s followed by a "1" ap- 
pears exactly once in the result (assuming cyclic boundaries). 
From the sequence's position, one easily infers t. For defi- 
niteness, assume t = 0. We then encounter a cluster state in 
the encoding \Ok) and |Wfc). The key point to notice is that, 
by Ref. |2lll . any two pure orthogonal multi -partite states can 
be deterministically distinguished using LOCC. In fact, em- 
ploying the construction of Ref. 112 111 , this can be done by 
an appropriate ordered sequence of adapted projective mea- 
surements 7Ti <S> ■ ■ ■ ® Ttk on the sites of each codeword, the 
effect corresponding exactly to an arbitrary given projective 
dichotomic measurement with Kraus operators IV'XV'I and 



IV* /(V* I = 1 — ("01 in the logical space. Hence, one can 
translate any single-site measurement on a cluster state into 
an LOCC protocol for the encoded cluster. This shows that 
|\&) is universal for deterministic MBC. At the same time, the 
von Neumann entropy S v n of any site of the initial resource 
is arbitrarily small for sufficiently large k: From the distri- 
bution of "0"s and "l"s in the codeword, one finds that the 
entropy for a measurement in the computational basis reads 
Sz = -ff&(3/(4fc + 2)), where is the binary entropy func- 
tion. Using the concavity of the entropy function, we have 
that LE < S v n < Sz- It follows that not even a Bell pair can 
be deterministically created between any two fixed systems. 

Outlook. - Until now, the only known scheme for MBC 
was the QCc and slight variations. Entire classes of states 
with physically reasonable properties (e.g. non-maximal lo- 
cal entanglement, long-ranged correlations, weakly entangled 
bonds) could not be dealt with. It is our hope that the frame- 
work presented opens up the possibility to adopt the compu- 
tational model to some extent to the specific physical systems 
at hand and no longer vice versa. For example, in linear optics 
computing, bonds are the easier to create the lower the entan- 
glement 12311 . Under those circumstances, there may well be a 
trade-off between the effort used to prepare a resource and its 
efficiency for MBC JH]. In turn, for cold atoms in optical lat- 
tices, exploiting cold collisions l24ll . configurations as in Fig. 
1 (a) could as feasibly be created as the cluster state, mak- 
ing use of a different interaction time for diagonal collisions. 
Other states may well be less fragile to finite temperature and 
decoherence effects. The presented tools open up a way for 
studies of quantitatively exploring such trade-offs in prepara- 
tion. 
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